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ABSTRACT 



Context. NGC 5253 is a nearby peculiar Blue Compact Dwarf (BCD) galaxy that, on account of its proximity, provides an ideal 
laboratory for detailed spatial study of starburst galaxies. An open issue not addressed so far is how the collisional and self-absorption 
effects on He i emission influence the determination of the He + abundance in 2D and what is the relation to the physical and chemical 
properties of the ionized gas. 

Aims. A 2D, imaging spectroscopy, study of the spatial behavior of collisional and radiative transfer effects in He + and their impact 
on the determination of He + abundance is presented for the first time in a starburst galaxy. 

Methods. The He i lines are analysed based on previously presented optical Integral Field Spectroscopy (IFS) data, obtained with 
FLAMES at the VLT and lower resolution gratings of the Giraffe spectrograph as well as with GMOS at Gemini and the R381 grating. 

O Results. Collisional effects reproduce the electron density structure. They are negligible (i.e. ~0. 1-0.6%) for transitions in the singlet 
cascade while relatively important for those in the triplet cascade. In particular, they can contribute up to 20% of the flux in the 
He 1/17065 line. Radiative transfer effects are important over an extended and circular area of ~30 pc in diameter centered at the 
Super Star Clusters. The singly ionized helium abundance, y + , has been mapped using extinction corrected fluxes of six Hei lines, 
realistic assumptions for T e , n e , and the stellar absorption equivalent width as well as the most recent emissivities. We found a mean(± 
standard deviation) of 10 3 >' + ~ 80.3(±2.7) over the mapped area. The relation between the excitation and the total helium abundance, 
y fo „ is consistent with no abundance gradient. Uncertainties in the derivation of helium abundances are dominated by the adopted 
5^ ' assumptions. We illustrated the difficulty of detecting a putative helium enrichment due to the presence of Wolf-Rayet stars in the 

main GH nR. Data are marginally consistent with an excess in the N/He ratio in the nitrogen enriched area of the order of both, the 
atmospheric N/He ratios in WR stars and the uncertainties estimated for the N/He ratios. We explored the influence of the kinematics 
in the evaluation of the Hei radiative transfer effects. Our data empirically support the use of the traditional assumption that motions 
in an extragalactic Hn region have a negligible effect in the estimation of the global optical depths. Individually, the broad kinematic 
component (associated with an outflow) is affected by radiative transfer effects in a much more significant way than the narrow one. 
^ . We found a relation between the amount of extra nitrogen and the upper limit of the contribution from radiative transfer effects that 

C^*l ' requires further investigation. We suggest the electron temperature as perhaps a common agent causing this relation. 
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1. Introduction the su pernebula in whic h gas is outflowing from the molecular 

__l ■ , : | , cloud ( Beck et afl. l2012h . Indeed, the whole central region of the 

■ At a distance of 3.8 Mpc ( gakai et al| | 2004f), NGC 5 253J galaxy is dominated by an intense burst of star formation in the 
> ; the Centaurus A / M 83 group d Karachentsev et al| |2007j), form of a large collect ion of compact young (-1-12 Myr) star 
£ . is one of the closest Blue Compact Dwarf (BCD) galaxies. dusters ( | H arris et al.l |2004). In contrast to this, the main 
l£j ; This galaxy is well known for presenting several peculiarities, body of NGC 5253 resembles mat of a dwarf elliptical galaxy 
, whose detailed study are closely connected to its proximity and and recentlyi ^ potentially massive ( > lfJ 5 Mq) and old (1 _ 2 
high surface brightness. For example, it contains a deeply em- Gyr) star dusters haye been found in the outskirts of the galaxy 
bedded very dense compact H n region at its nucleus (here- (Harbeck e JMM^I- Finally, NGC 5253 is best-known for be- 
after, "the supemebukfX detected in the radio at 1 .3 cm and ing Qne of ^ few examples (and the closest) of a galaxy present . 
2 cm ([Turner et aU 2000|) that host two very massive Super ing a confirmed local excess in nitrogen (see e g . Walsh & Rov> 
Star Clusters (SSCs. lAlonso-Herrero et al.l 12004) and is embed- 1939) 
ded in a larger (i.e. ~100 pcx80 pc) Giant Hn Region (here- 
after, the central GH11R). Recently, mid-infrared observations We "re carrying out a detailed study of this galaxy 
showed how its kinematics is compatible with a model for usln 8 Inte § ral Fleld Spectroscopy (IFS). The results ob- 
tained so far have f urther highlighted its peculiar nature. In 
Send offprint requests to: A. Monreal-Ibero |Monreal-Ibero_etalJ 420101 hereafter Paper I), we found that 
* Based on observations collected at the European Organisation the emission line profiles were complex and consistent with 
for Astronomical Research in the Southern Hemisphere, Chile (ESO an scenario wher e the two SSCs produce an outflow (see also 
Programme 078.B-0043 and 383.B-0043) and at the Gemini South iBeck et all |2012|) . Also, we delimited very precisely the area 
Telescope (Programme GS-2008A-Q-25). polluted with extra nitrogen. Moreover, we detected nebular 
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He H/14686 in several locations, some associated with WN-type 
Wolf-Rayet (WR) stars (as traced by the blue bump at around 
4680A) and some not, but not nec essarily coincident w i th the 
area exhibiting extra nitrogen. In Monrea FIbero et aD ([2012, 
hereafter Paper II), we studied the 2D distribution of the main 
physical (electron temperature and density, degree of ioniza- 
tion) and chemical properties (metallicity and relative abun- 
dances of several light elements) of the ionized gas. A new 
area of enhanced nitrogen abundance at ~130 pc from the main 
area of enhancement and not reported so far was found. In 
IWestmoquette et ail d2013l hereafter Paper III) several locations 
showing emission characteristic of WC-type WR stars (via the 
red bump at around 5810A) were identified. The fact that WR 
stars are spread over ~350 pc gives an idea of the area over which 
the recent starburst has occurred. The chemical analysis was ex- 
tended with the finding that, with the exception of the aforemen- 
tioned localised N excess, the O/H and N/H distributions are 
flat within the whole central 250 pc. 

An issue not addressed in detail so far in NGC 5253 is 
the 2D determination and distribution of the He + abundance. 
In a cosmological context, this is particularly relevant since 
the joint determination of metallicity (as traced by the O/H 
abundance) and 4 He abundance (}fl) for extragalactic Hn re- 
gions and star-forming galaxies at low-metallicity was proposed 
as a means to estimate the primordial helium abund ance, Yp 
dPeimbert & Torres-Peimbertl 119761 iPagel et all 119861) serving 
as a test-bench for the standard hot big band model of nucle- 
osynthesis. However, the density of baryonic matter depends 
weakly on Yp. Therefore, to put useful constrains on Yp, A He 
abundance of individual objects has to be determined with ac- 
curacies <1%. Nowadays, emission line flux data of this qual- 
ity can be achieved and, indeed, the astronomical commu- 
nity is actively work ing on gett i ng and improving the estima- 
tion of Yp (see e.g. lAver et al.L l2010t llzotov & ThuanT, 12 010: 
Fukugita & Kawasaki, 2006; Peimbert et all 120071: llzotov et al.L 
2007tloiive & SkillmanLl200lHPeimbert et all l2002l for recent 
estimations by the different groups). However, He abundance de- 
terminations are influenced by several effects and systematic er- 
rors which are not, in pri nciple, straightforward t o quantify and 
correct (see for example lOlive & Skillmanl l200lb . Specifically, 
the intensity of Hei emission lines may intrinsically deviate 
from the recombination values due to collisional and radiative 
transfer effects. Moreover, the emitted spectrum also depends 
on the physical conditions of the ionized gas (e.g. temperature, 
density and ionization structure). Also, on top of these specific 
properties of the Hn region, extinction by dust and a possi- 
ble underlying stellar absorption component can also affect the 
observed spectrum. All these effects contribute to the uncer- 
tainties associated with the estimation of ionized helium abun- 
dance (y + = He + /H + ). A final extra source of uncertainty is 
associated with the estimation of the amount of existing neu- 
tral helium (i.e. the estimation of the ionization correction fac- 
tor, icf(He)) and the calculation of the total helium abundance, 
y = y to , = icf(He) x y + . 

Tentative values for y + were presented in Paper I based on 
the He 1/16678 line, which is almost insensitive to collisional and 
self-absorption effects. However, 2D distributions of all the rel- 
evant physical properties for the ionized gas were not available 
at that time. Moreover, other He i lines, in particular He 1/17065, 
can indeed be affected by collisional and self-absorption effects, 



1 Here, we use Y for the helium mass fraction and y for the number 
density of helium relative to hydrogen. Assuming Z = 20(0/ H), they 
are related as K = Mk*Mm , 

l+4y 



specially in conditions of relatively high electron temperature 
(T e ) and density (n e ), as it is the case in the main GHuR of 
NGC 5253. Supporting this expectation, long-slit measurements 
predict too high He + abundances from the He 1/17065 when only 
recombination effects are taken into account dSidolil [2010h and 
a non-negligible optical depth r(3889) when the He + abundance 
is derived using several em ission lines in a consistent manner 
(Lop ez-Sanchez et alll2007l) . 

At present, a complete 2D characterization of the physical 
properties of the ionized gas in the main GH iiR of NGC 5253 
is available. Therefore, we are in an optimal situation for both 
mapping the collisional and radiative transfer effects in 2D, and 
re-visiting the derivation of the He + abundance map taking into 
account many Hei lines. This will be the main purpose of this 
work . Given t hat the metallicity of the object (12 + log(0/H) = 
8.26, Paper II) is only moderately low, our focus will not to be in 
achieving the < 1 % accuracy required in the determination of the 
primordial He abundance, but to explore the effects of a param- 
eter not taken into account so far: namely spatial resolution. In 
addition, we will explore the 2D relationship between the prop- 
erties of the ionized gas derived so far and the He i collisional and 
self-absorption effects. To our knowledge this work constitutes 
the first attempt to study in 2D the collisional and self-absorption 
effects in He i in any extra-galactic object. Moreover, irrespective 
of the spatial resolution, due to the characteristics inherent in IFS 
data, there is the guarantee that the set of ~ 100 spectra utilized 
in this work have been processed in a homogeneous manner all 
the way from the observations (i.e. a given observable was taken 
with the same observing conditions for the whole set of data) to 
the final y + and y tot derivations. 

The paper is structured as follows: Sec. 2 describes the char- 
acteristics of the data utilized in our analysis; Sec. 3 contains an 
evaluation of the He i collisional and self-absorption effects in 
2D as well as the derivation of the y + and y maps. Sec. 4 dis- 
cusses the relation between radiative transfer effects and other 
quantities (e.g. kinematics of the gas, relative N/O abundance). 
Our main conclusions are summarized in Sec. 5. 



2. The data 

We will focus our study on the area associated with the main 
GHuR in NG C 5253 (see Fig. Q}. This is a portion of the full 
area studied in Paper I and the location where: i) the gas presents 
relatively high electron temperature and density, and therefore, 
important collisional and radiative transfer effects are expected; 
and ii) several helium lines can be detected over a relatively 
large area with sufficient quality, and therefore a 2D analysis 
based on the information in individual spaxels is feasible. The 
utilized data were collected during several observing runs using 
the FLAMES-Argus and the GMOS Integral Field Units (IFUs) 
and together cover the whole optical spectral range. In the fol- 
lowing, we briefly describe the basic instrumental characteris- 
tics of each set of data and compile the information that was 
extracted. 

2.1. FLAMES-Argus data 

Data were obtained with FLAMES dPasquini et all l2002h at 
VLT UT2 in Paranal. We used the Argus IFU with the sampling 
of 0.52"/lens, covering a field of view (f.o.v.) of 1 175x7'.'3, and 
four low resolution gratings (LR1, LR2, LR3, and LR6). All to- 
gether, they offer a spectral coverage of 3 610 - 5 070 A plus 
6440-7 180 A at a dispersion of 0.2 Apix~'. Details of the ob- 
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Fig. 1. Left: False colour image in filters F658N (Ha, cyan 
channel), F55QM (V, yellow channel), and F814W (/, magenta 
channel) for the central part of NGC 5253 using the HST-ACS 
images (programme 10608, P.I.: Vacca). The area studied here is 
marked with a black rectangle. Right: Ionized gas distribution as 
traced by the extinction corrected H/3 map derived from a por- 
tion of the original FLAMES data. The position of the two main 
peaks of continuum emission are marked with crosses. The map 
is presented in logarithmic scale in order to emphasize the rele- 
vant morphological features and cover a range of 2.0 dex. Flux 
units are arbitrary. Note the existence of three dead spaxels at 
~ [5'.'0, -l'.'O] as well as absence of signal in the spaxels at the 
two left corners of the field of view. 

servations, data reduction, and cube processing, as well as maps 
for most of the physical and chemic al proper ties utilized in this 
work, have already been presented in lPaper il and lPaper III Maps 
for the helium lines were derived by independently fitting in 
each spaxel the He i line profiles with a single Gaussian function 
with MPFITEXPR ( Mar kwardtL [20091) . For the particular case of 
Hei/13889, which is blended with H8 at our spectral resolution, 
we created an extinction corrected map by subtracting from the 
H8+Heivi3889 map, a map of 0.659x the H7 map. Thi s H8 line 
intensity is that predicted bv lStorev & Hummer! (1 1995b for Case 
B, T e = 10 4 K and n e = 100 cm' 3 . The final set of FLAMES- 
Argus maps utilized here are: 

i) an extinction map derived from the Ha/H/3 line ratio; 

ii) a map for the H/3 equivalent width (EW(H/3)); 

iii) a map of electron temperature T e as derived from the 
[Oiii]^4959,5007/ [Ora]i4363 line ratio: r e ([Ora]). In 
those spaxels where no determination of r g ([O ml) was 
available, we assumed r e ([Om]) = 10500 K (see iPaperll 
for typical r e ([Om]) values outside the main GH nR); 

iv) a map for the electron density {n e ) as derived from the 
[Sii]/16717/[Sii]/}6731 line ratio; 

v) maps for the + /H + , O/H and S + /H + abundances, as de- 
rived from collisional lines using the direct method, to esti- 
mate the icf(He); 

vi) maps for different tracers of the excitation degree (i.e. 
[Sn],M6717,6731/Ha and [Om]/I5007/Hy3 line ratio) to be 
used in the estimation of the icf(He) at those locations 
where no measure of the + /H + , O/H and S + /H + abun- 
dances is available and to explore the dependence of y + and 
y m on the excitation; 

vii) a map of the relative abundance of nitrogen, N/O, as derived 
from collisional lines using the direct method; 

viii) maps for the /T3889, ^4471, /14922, ^6678, Z7065 Hei 
equivalent widths and extinction corrected line fluxes using 
our ex tinction map and the extinction curve of Fluk s et al.l 
(11994 . The extinction corrected line flux maps normalized 
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Fig. 2. Extinction corrected flux maps normalized to H/3 for the 
lines utilized in this work. The position of the two main peaks 
of continuum emission are marked with crosses in this and all 
subsequent maps. We marked in white the areas where a given 
line was not observed. Specifically, for the He 1/16678, He 1/14922 
and He 1/17065, they correspond to dead fibers. For He 1/15876, 
He 1/14471 and He 1/13889, these areas were not covered by the 
corresponding GMOS or FLAMES-Argus field of view. 



to H/3 are presented in Fig. [2] Note that in order to mini- 
mize uncertainties associated to aperture matching, absolute 
flux calibration and extinction, lines were measured rela- 
tive to a bright Balmer line observed simultaneously with a 
given helium line. Then, we a ssumed the theoretical B almer 
line intensities obtained from I Storey & Hummer! (1 19951) for 
Case B, T e = 10 4 K and n e = 100 cnT 3 . Additional Hei 
lines were covered by the FLAMES set-up but not used 
here. Specifically, He 1/15016 and He 1/14713 emission lines 
were detected overmost of the FLAMES f.o.v. He 1/15016 is 
relatively close to the much brighter (i.e. ~ 150-400 times) 
[Om]/i5007 line and, in the main GHiiR, the wings of the 
[Ora] line profile prevented us from measuring a reliable 
line flux. Regarding the relatively weak /14713 line, at the 
spectral resolution of these data, this line is strongly blended 
with [Ariv]/i471 1 and the uncertainties associated with the 
deblending of these lines in the GH nR are relatively large 
due to th e existence of seve ral distinct kinematic compo- 
nents (see Paperll: [Paper IIll) : 
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ix) maps with the extinction corrected line fluxes in /16678, 
/i7065 He i and Ha, for the different kinematic components 
presented in lPaperll 

2.2. GMOS data 

The Gemini-South Multi-Object Spectrograph (GMOS) 
data were taken us i ng th e one-slit mode of its IFU 
dAllington-Smifh et all |2002|) . In this mode, the IFU cov- 
ers a f.o.v. of 5'.'0 x 3'.'5 sampled by 500 contiguous hexagonal 
lenslets of 072 diameter. The utilized grating (R381) gives 
a spectral coverage of 4 750 - 6 850 A at a dispersion of 
0.34 A pix ' 1 , thus c omplementing the FLAMES-Argus data. 
We refer to Pa per Ilj for further details on the observations and 
data reduction. The product of the reduction is a datacube per 
pointing with a uniformly sampled grid of 07 1 . In this work, we 
utilized the two pointings (out of four) that mapped the central 
GHnR. As we did with the FLAMES's data, in each spaxel 
all the lines of interest were independently fit with a single 
Gaussian function with MPFITEXPR. The final set of GMOS 
maps utilized here are: 

i) Maps for Ha and H/3 fluxes. These images were utilized to 
check the consistency between the FLAMES and GMOS 
data, both in terms of observed structure and derived ex- 
tinction map, to estimate the offset and rotation that was 
necessary to be applied to the GMOS data, and to correct 
for extinction the He 1/15876 map; 

ii) A map for He 1/15876 flux. This is the strongest He i line and 
therefore, one of the key observables for the present study; 

iii) Equivalent widths and extinction corrected line flux maps 
of He 1/15876 flux normalized to H/3. These were derived 
from the maps previously mentioned and were rotated and 
reformatted to match the FLAMES data using the drizzle 
task of the Space Telescope Science Data Analysis System 
(STSDAS) package of IRAlfl They were the only GMOS 
maps utilized jointly with the FLAMES's maps. The flux 
map is shown in Fig. [2] together with the other He i line flux 
maps. 

3. Results 

3.1. Collisional effects as traced by the theoretical C/R ratio 

Collisional excitation in hydrogen can be important in regions 
of very low m etallicities (Z<1 / 6 Z ) due to their relatively high 
temperatures (iLuridianai I2009I) . However, for the typical tem- 
peratures and densities found in Hn regions in general, and 
in the main Giant Hn Region (GHuR) of NGC 5253 in par- 
ticular, the collisional excitation in hydrogen is negligible in 
comp arison with recombination (see e.g. Fig. 9 of lAver et all 
2010). Thi s is not the case for he lium. The Hei 2 3 5 level (see 
Fig. 1 in iBeniamin et all [1 999, for a representation of the 
Grotian diagram for He i singlet and triplet ladders) is highly 
metastable, and collisional tr ansition from it can be important 
(lOsterbrock & Ferlandl 120061) . Specifically, at relatively high 
densities this level can be depopulated via collisional transitions 
to the 2 3 P°, 2 1 P° and 2 l S and, to a lesser extent, to higher sin- 
glets and triplets (mainly 3 3 P°). The effect on the observed Hei 

2 The Image Reduction and Analysis Facility IRAF is distributed by 
the National Optical Astronomy Observatories which is operated by the 
association of Universities for Research in Astronomy, Inc. under coop- 
erative agreement with the National Science Foundation. 



emission lines is always an increase in the observed flux. The 
relative importance of the collisional effects on a given emission 
line is characterized by the C/R factor, i.e. the ratio of the col- 
lisional component to that arising from recombination which is 
given by: 

£ = n^gkeff 

R n He a eff 

where «2 3 s> an d n He are tne densities of the 2 3 S state and 
He + respectively, a e ff is the effective recombination coefficient 
for the line, and k e ff is the effective collisional rate coefficient. 




Ax (arcsec) Ax (arcsec) 

Fig. 3. Collisional effects as traced by the C/R ratio for the He i 
lines utilized in this work. Note that a common scale was used 
for all the lines in order to facilitate comparison of the rela- 
tive effects between the different lines. Also, a logarithmic color 
stretch is used to emphasize the variations within the region for 
a given line. 

Here, we will estimate the 2D contribut ion to the co l lisiona l 
component using the relations derived by Port er et aD (120071) . 
where theoretical C/R factors are calculated as functions of n e 
and T e , and assuming that the density from [S n] traces well that 
from helium. Plasma temperatures as traced by helium lines can 
typic ally be ~50% of th ose traced by oxygen in planetary nebu- 
lae (|Zhang et afl, l2005l) . For H n the situation is not so clear. For 
the specific case of NGC 5253, values between 82% and 96% 
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Fig. 4. Maps for the ratio between the C/R factors derived for 
two assumptions of the T e (He i): Q(A) = Note that ar- 

eas with constant values for which a r f ([Om])=10500 K was 
assumed have not been included in the comparison. A logarith- 
mic color stretch is employed to emphasize the variations within 
the region for a given line. 



have been reported in specific apertures (Lope z-Sanchez et al.L 
l2007l) . We made two assumptions for the electron temperature: 



Secondly, lines corresponding to transitions in the singlet 
cascade have a negligible contribution from collisional effects 
(e.g. C/R factor for A4922 varies between ~ 0.001 - 0.006) 
while for those lines in the triplet cascade (A1Q65, /15876, /14471, 
and /13889), the contributions from collisional effects can be im- 
portant. In particular, the C/R factor for /i7065 ranges between 
~ 0.02 and ~ 0.22, meaning it reaches ~20% in the nucleus of 
the galaxy. 

Thirdly, the assumed temperature has some influence in the 
estimation of the collisional effects. In our particular case the 
assumed temperature in Case 2 was only ~ 15% smaller than 
that in Case 1 . However, this implies a smaller contribution of 
the collisional effects by -25-30 % for A1Q65 and -30-35 % for 
/13889 (the two lines most affected by collisional effects) and up 
to ~ 50% for the other lines under study. Interestingly, areas of 
lower temperature are more sensitive to the assumption on T e . 
It is important to note that the uncertainties associated to the 
errors due to the measurement of the line fluxes invol ved in th e 
determination of r e ([Om]) were typically <1 000 K (PaperB. 
This is smaller than the difference in the assumed temperature 
between the two reasonable assumptions, Case 1 and 2, which 
range between ~1 500 K at the peak of emission to ~1 300 K in 
the areas of lowest surface brightness. This implies that, at this 
level of data quality, more than errors associated to line fluxes, it 
is systematic errors associated with assumptions in density and 
temperature that dominate the uncertainties in the evaluation of 
the collisional effects. 



3.2. Radiative transfer effects as traced by r(3889) and 
derivation of He + abundance 

Single ionized helium abundance, y + , can be calculated as fol- 
lows: 



y + W) = 



F(A) E(H/3) EW( ew { ")' W 1 
F{Hp) E{A) W^Mim/AA) 

tW(A) 



(2) 



where F(A)/F(H/3) is the extinction corrected flux for the 
He i line scaled to H/3, E(H/3) and £(Hei) are the theoretical 
emissivities; f T (A) is a factor that takes into account radiative 
transfer effects. The remaining term in Eq. [2] with a(A) is the 
equivalent width in absorption, which takes into account the ef- 
fect of the underlying stellar population. In the following, we 
will describe how each of these terms were evaluated. 



- Case 1 : The temperature from oxygen traces well that from 
helium: T e (He i) = T e ([0 m]); 

- Case 2: The temperature from helium is proportional to the 
temperature from oxygen, with the constant of proportional- 
ity estimated as the mean of the ratios be t ween b oth temper- 
atures provided by lLopez-Sanchez et aH (120071) : T e (He i) = 
0.87 T e ([Q m]). 



3.2.1 . Calculation of the emissivities 

For H/3 emissiv it ies, we utilized the function from 
IStorev & Hummed d 1995b which is, in units of 10 erg 



E(Hff) = 4irj H p/n e n H + 



l.31tf m exp(-0.104/f e ) 



(3) 



Maps for the C/R factors using the first assumption are pre- 
sented in Fig. [3] while the ratio between the two estimations is 
presented in Fig. [4] This is the first time that collisional effects 
for a set of helium lines have been mapped in an extra-galactic 
source. Several results can be extracted from these maps. 

Firstly, all the C/R maps display the same structure. That 
is: higher ratios at the peak of emission for the ionized gas and 
towards the northwest half of the GH nR and a decrease of the 
collisional contribution outwards. T his repro duces the observed 
density structure (see e. g. Fig. 6 in iPaper lb . 



with t e = T e /1Q 4 . For the He i lines, we utilize d those emis- 
sivities origina lly provid e d bv [ Porter et al.l (1201 2b and recently 
corrected bv lPorter et al.l (12013b . These are the most recent He i 
emissivities and have the collisional effects included. Therefore, 
there was no need to include an extra term in Eq. |2]to take into 
account collisional effects. They are tabulated for discrete values 
of n e and T e . However, the n e in the main GH nR varies between 
<100 cirT 3 and ~ 660 cirT 3 and T e ranges between ~ 9 000 and 
~ 12 000 K. Therefore, to evaluate the emissivities at each in- 
dividual spaxel, we fitted the values provided for n e = 100 and 
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1 000 crrT 3 and T e ranging from 5 000 K to 25 000 K to functions 
with the same parametrization as E(H/3), a^ e exp(c/f f ), and then 
interpolated on a logarithmic scale as follows: 



E(A, log n e ) = E(A, 2) + (E(A, 3) - E(A, 2))(log n e - 2) 



(4) 



The coefficients and standard deviations of the fits are com- 
piled in Table Q] while Fig. [5] shows a comparison between the 
fitted function and the discrete values provided by iPorter et"al] 
(2012) for the range of densities and temperatures covered in the 
GHnR. A comparison between this figure and Fig. [3] shows the 
correspondence between the degree of dependence on the den- 
sity and the contribution of the collisional effects. 




Fig. 5. Fitted functions to the IPorter et al.l (1201 2|) emissivities. 
The black lines are the functions for n e — 100 crrT 3 (continuous 
line) and for n e = 1000 cirT 3 (long dashed line). The intemedi- 
ate short dashed lines represent the interpolated emissivities for 
log(n e ) =2.2,2.4,2.6,2.8. 



3.2.2. Correction for underlying stellar population 

To take into account the effect of the underlying stellar pop- 
ulation it is necessary to estimate the equivalent width both 
in emission and absorption for H/3 and the helium lines. The 
equivalent width of H/3 in emission (not shown) ranges typ- 
ically from -240 A at the peak of emission for the ionized 
gas to ~65 A in the most outer regions of the area sam- 
pled here. We assumed a component in absorption of 2 A, 
which is adequate for very y oung starbursts, as the one at 
the nucleus of NGC 5253 ( e.g. iGonzalez Delgado et all 12005b 
lAlonso-Herrero et all 1201 0b . This implies a correction for H/3 
in absorption from a negligible value (i.e. <1%) at the peak of 
emission to ~4% in the outer parts of the covered area. 

The correction due to the underlying stellar population is not 
st raightforw ard to estimate for the helium lines. As explained 
in Pa per III for the data obtained with FLAMES Giraffe and the 



Table 1. Coefficients for the equations fitted to the IPorter et al.l 
(120121) Table 2. 



Line 


a 


b 


c 


Std Dev (%) 


n e = 100 cnT J 


3889 


1.453 


-0.724 


-0.036 


0.010 


4471 


0.679 


-1.078 


-0.105 


0.003 


4922 


0.184 


-1.091 


-0.105 


0.001 


5876 


1.680 


-1.061 


+0.004 


0.042 


6678 


0.473 


-1.065 


+0.013 


0.015 


7065 


0.269 


-0.367 


+0.100 


0.007 


n e = 1000 cm^' 


3889 


1.198 


-0.336 


+0.207 


0.032 


4471 


0.507 


-0.694 


+0.202 


0.048 


4922 


0.131 


-0.642 


+0.252 


0.016 


5876 


0.895 


-0.195 


+0.671 


0.255 
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0.237 


-0.138 


+0.738 


0.094 
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0.320 


+0.171 


+0.163 
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Fig. 6. Maps of the correction factor due to a component in ab- 
sorption in the He i lines. Different scales were utilized in the dif- 
ferent maps to emphasize the structure. The corresponding map 
for the H/3 correction factor (not shown) is similar to those for 
,16678 and ,17065. 



LR1 or LR2 gratings, the contribution of the stellar population 
was separated from that of the gas using the STARLIGHT code 
dCid Fernandes et afl 120051 [2009) to match the stellar contin- 
uum. Therefore, we can assume a// e ,/(,14471) = fljy c /(^3889) = 
0. For the other lines, a common approach assumes identical 
stellar equivalent widths for all the helium lines. However, Hei 
lines are produced by the bluest stars. Therefore, in the con- 
text of a galaxy suffering a burst of star formation on top of 
an older population, redder lines should have smaller equiva- 
lent widths in absorption, since older stars contribute more to the 
stellar continuum. Specifically, typical estimations of the equiv- 
alent width in absor ption for redder lin es would be ~ 40 - 80% 
that of aHei(/l4471) (lAver et all 120101) . This is one of the lines 
observed with the LR2+FLAMES configuration for which neb- 
ular and stellar information have been disentangled. Therefore 
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a// el (/14471) could be measured from our emission line free 
cube. Typical values were extremely low (i.e. ~ 0.11 + 0.03 A) 
and without any obvious variation following the structure of the 
GHuR or the location of the Super Star Clusters. Taking this as 
reference, a nd the relative valu es between the different a//ei(^)'s 
reported by I Aver et al. (2010 |, which were in turn derived us- 
ing the models prese nted by Gonzalez Delgado et alJ (120051) and 
iMartrns et ail (120051) . we assumed a HeI (X) = 0.09, 0.08, 0.06 and 
0.05 A for ,14922, ,15876, ,16678, ,17065, respectively. Equivalent 
widths for the ,14922, ,15876, ,16678, ,17065 Hei emission lines 
were ~ 0.3 - 4 A, ~ 20 - 60 A, ~ 3 - 20 A, and ~ 3 - 35 A, 
respectively. With these values, the largest correction for absorp- 
tion was for He 1/14922, with typical values between 3% and 6%. 
However, corrections in the most external spaxels could reach up 
to ~25%. For the other lines the correction was more moderate, 
with values between ~1% and <8% for He 1,15876 and compa- 
rable to the correction in H/3 for the He 1/16678 and He 1/17065 
lines. This is illustrated in Fig. [6] 

3.2.3. Estimation of radiative transfer effects 

Radiative transfer effects can be important in recombination ra- 
diation. Given the metastable character of the 2 3 5 level, un- 
der certain conditions, the optical depths in lower 2 3 S - n 3 P° 
lines imply non-negligible eff ects on the emission line strengths 
(Osterbrock & Ferland, 2006). Specifically ,110830 photons are 
only scattered, but absorbed photons, corresponding to transi- 
tions to higher levels, can be converted to several photons of 
lower energy. The most prominent example is ,13889 photons 
that can be converted to ,14.3 fim 3 3 S - 3 3 P°, plus ,17065 2 3 S - 
3 3 P°, plus ,110 830 2 3 S -2 3 P°. The net effect is that lines associ- 
ated with transitions from the 2 3 P level upwards (e.g. ,13889) are 
weakened by self-absorption, while lines associated with several 
transitions from higher levels (e.g. ,17065) are strengthened by 
resonance fluorescence. Contrary to collisional effects, radiative 
transfer effects do not affect photons in the singlet cascade. The 
relative importance of radiative transfer effects is quantified by a 
correction factor, f T (A), for each line which is a function of the 
optical depth at ,13889, t(3889). Here, we parametrized these 
factors by fitting the values provided by Robbins (1968) to a 
non-expanding nebula to the functional form / T (/l) w= o = 1 + aT b . 
The corresponding functions for the lines utilized in this work 
are: 



/r(7065) 6 
A(6678) 6 
/r(5876)„ 
/r(4922) 6 
/r(4471)„ 
/t(3889) 6 



1 + 0.741T 0341 
1 

1 + 0.0126T 0496 
1 

1 + 0.0022-r - 728 



1 -0.26 It 



.0.305 



(5) 
(6) 
(7) 
(8) 
(9) 
(10) 



Note that no values were pr ovided for He 1,14922 in the 
original work of iRobbinsI (11968). However, as it happens for 
He 1/16678, this is a singlet line, and therefore / T (4922) = 1 can 
be assumed. 

Typically, optical depth (r(3889)) and He i abundance (y + ) 
(and other parameters) are determined simultaneously by mini- 
mizing^ 2 , defined as the difference between each helium line's 
abundance (weighted according to a reasonable criterion like the 
line flux) and the average. In this methodology it is implicit 



that all the lines trace the same location in the nebula/galaxy. 
However, the area under study in this work suffers from heavy 
extinction (see Fig. 3 in iPaperfl) . Therefore, a priori, it is not 
possible to assume that all the lines equally penetrate the nebula 
interior and that bluer and redder lines trace zones with the same 
y + . Because of that, we grouped the lines in two sets according 
to their wavelengths, called hereafter the blue (,13889, ,14471, 
,14922) and red (,15876, ,16678, ,17065) sets, which will be ana- 
lyzed independently. Each set is made out of: i) a line from the 
singlet cascade, and therefore not affected by radiative transfer 
effects; ii) a line from the triplet cascade, thus highly sensitive 
to radiative transfer effects; iii) a line from the triplet cascade 
mildly sensitive to radiative transfer effects. 
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Fig. 7. Maps for y + derived from lines of the singlet cascade for 
the blue (left) and red (right) sets. A common scale was used in 
both maps to facilitate comparison between them. 

A first estimation of the y + abundance structure as traced by 
the red and blue lines can be obtained from the singlet lines, 
since they are not affected by radiative transfer effects. This is 
shown in Fig. [7] The mean (± standard deviation) are 80.7(+5.1) 
and 76.8(+1.8) for the ,14922 and ,16678 lines averaged over the 
mapped area. These values indicate that even if the red and blue 
lines are not tracing exactly the same gas columns, at least they 
sample areas with a the same y + within ~5%. 

A comparison of the initial y + maps with those obtained from 
a line highly sensitive to radiative transfer effects in each set 
(,13889 and ,17065), allowed us to determine the respective f T (A) 
map for each set, which can, in turn, be converted into r(3889) 
maps. These are shown in Fig. [8] The optical depth has the same 
structure in both maps, with the peak at the main super star clus- 
ters). The shape of the area presenting t(3889) > is circular 
but resolved. With a FWHM= 176~30 pc, this is larger than the 
seeing (~079). This is consistent with a picture where large op- 
tical depths are not restricted to the deepest core of the galaxy, 
associated with the supernebula, but extend over a larger region. 
Also, the fact that the optical depths derived from the red set of 
lines are larger than those derived from the blue set is in har- 
mony with a picture where deeper layers in the nebula, which 



we interpreted as denser and hotter (Monreal-Ibe ro et al.i 1201 
see also Beck et al. 2012), suffer from larger radiative transfer 
effects and suggests a link between He i optical depth and dust 
optical depth. 

3.2.4. Final derivation of He + abundance 

Finally for both sets, the information of the mildly sensitive line 
was added and abundances from each line were recalculated 
using their corresponding t(3889). The final abundance maps 
for each set were made using a weighted average. We used the 
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Ax (arcsec) Ax (arcsec) 

Fig. 8. Maps for r(3889) as derived from the /14922 and ,13889 
(left) and the /16678 and /T7065 (right) emission lines. 



mean fluxes of the He i in the utilized area to determine the re- 
spective weights. These were 3:1:1 for the ^5878:^16678:^7065 
and 25:5:1 for ^3889:^447 1:^4922. The mean (+ standard de- 
viation) for the red and blue sets for y + are 79.3(+2.5) and 
82.0(+3.8), i.e., results from the red and the blue sets agree to 
within ~3%. Since differences between the maps derived from 
the blue and red sets are consistent with tracing the same y + per 
spaxel, we derived a final map using the information of all the 
lines by averaging these two maps with a blue:red weight of 
0.8:1.0, the mean ratio between the brightest lines in the blue 
and red sets (i.e. -13889 and /15876). This is shown in Fig. [9] 
The mean (+ standard deviation) is 80.3(±2.7). This value com- 
pares wel l with other value s of y + reported in the literature for 
this area dKobulnickv et all 1 19971 lLopez-Sanchez et all 120071: 
ISidoil2010ir ~ 
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Fig. 9. Final map fory + derived from a weighted average of those 
derived for the lines in the blue and red sets. Note we used a dif- 
ferent scale than in Fig.Qto increase the contrast and emphasize 
the structure following the excitation structure. Instead, the scale 
is common with that in the right column of Fig.[T2]to make eas- 
ier the comparison between total and ionic helium abundances. 



Finally, it is worth mentioning that throughout all the deriva- 
tion of y + , we have assumed that He i singlets are formed under 
Case B conditions (i.e. lines are formed in the limit of infinite 
Lyman line optical depth and there is no optical depth in tran- 
sitions arising from excited states). This is the standard frame- 
work for the derivation of helium abundances. However, under 
certain conditions, helium Lya /1584 (and higher Lyman transi- 
tions) may be attenuated due to two effects: abso rption of t he he- 
lium Lyq line by dust and hydrogen absorption dFerlandi fl980; 
Shield! fT993h . If these effects were important, the Case B as- 
sumption would no longer be applicable. The fact that we have 
recovered similar helium abundances from the two sets of lines 
supports the Case B assumption. 

3.3. Determination of He abundance 

The favorite scenario to explain the extra nitrogen found within 
the area studied here is that this has been produced by Wolf- 
Rayet stars and presumedly those at the two Super Star Clusters. 
In this scenario, as it happens with the nitrogen, an overabun- 
dance of helium is also expected to be observed. Specifically, 
typical stellar atmospheric N/He ratios would be ~ 3 x 10~ 3 and 
~ 5 x 10~ 4 for WN and WC stars, respectively dSmith & W illis. 
[1981 . This could be use as a reference for the expected N/He 
ratios of the material newly incorporated into the warm ISM. Are 
our data supporting an enrichment in helium abundance compat- 
ible with these ratios and, therefore, supporting the Wolf-Rayet 
stars hypothesis? 




i i i i i i i i i i i i i 

0.6 0.7 0.8 0.9 
log([OIII]A5007/HP) 



Fig. 10. He + abundance vs. [Om]/l5OO7/H0. The first-degree 
polynomial fit is shown with a black line. The l-cr and 3-<x levels 
are marked with thick and thin red lines respectively. Mean and 
standard deviation of each 0.05 dex bin in log([Om]/15007/Hj8) 
are shown with blue diamonds and error bars respectively. The 
green horizontal dashed line shows the mean value. Data points 
corresponding to spaxels with log(N/0) > -1.3 have been 
marked with orange circles. 



To estimate the emissivities, the calculation of these abun- 
dances was made with assumed n e ([S n]) and r e (Hei) = T e (Hi) 
= r e ([Oin]). As for the mapping of the collisional effects, in or- 
der to evaluated the influence of the selection of T e , the deriva- 
tion presented here was repeated assuming T e (Hei) = T e (Hi) 
= 0.87 T e ([Om]) and finding a mean (+ standard deviation) of 
79.4(+3.0). This implies that the precise selection of T e has a 
small effect (~1%) in the determination of y + , in comparison to 
other factors like e.g. the correction for absorption or radiative 
transfer effects, as long as one keeps this selection within rea- 
sonable values. 



As can be clearly seen in Fig. [9] the y + map presents some 
structure that follows the excitation, even if the standard devia- 
tion fory + is small (~4% of the mean value). This is better s een 
in Fig. [10] which is an updated version of Fig. 14 in iPaper II but 
restricted to the area studied here. The figure presents y + for each 
individual spaxel versus a tracer of the excitation, in this case 
[Om]/I5007/H/?. For each bin of 0.05 dex, we overplotted the 
mean and standard deviation of y + values with blue diamonds 
and error bars and fitted all our data points with a first-degree 
polynomial. Considering that if we find < y + > in the bin of 
highest [Oiii]/15007/Hj8 to be larger than < y + > +cr(y + ) in the 
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bin of lowest [O in] xl5 007 /H/S indicates the presence of a gradi- 
ent, then this plot is consistent with a positive gradient in y + . 
Moreover, spaxels with an extra amount of nitrogen have, on 
average higher y + . Is this increase due only to variations in the 
ionization structure or are we witnessing an enrichment in the 
helium abundance? To answer this question, one needs to esti- 
mate and correct for the unseen neutral helium, removing in this 
way, the dependence on the excitation. Thi s is not straight for- 
ward as examplified by previous studies (e.g.lViegas et all [2000; 
ISauer & Jedamziktl2002trGruenwald et al.ll2002l) . 
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Fig. 11. Maps for the ionic abundances utilized in the derivation 
of the ionization correction factors. Left: + /H + . Right: S + /H + 

Within the spirit of keeping the analysis simple, here we 
compare the results using thre e approaches. The first two 



make use the icf 's proposed by Kunth & Sargent! d 1983b and 
iPeimbert & Torres- Peimbert (1 19771) : 



icf(He) K&S83 = (l-0.25O + /Or 1 
icf(He)p &TP 77 = (1 - 0.35 + /0 - 0.65 S + /Sy 1 



(11) 
(12) 



For the t hird one, we assumed a p riori a functional form sim- 
ilar to that of lKunth & Sargentl d 1983b and iteratively determined 
the required coefficients to obtain a relation consistent with no 
gradient of helium abundance. This would represent the case of 
largest reasonable correction for unseen helium: 



icf(He) M i3 = (1 -0.46O + /Or 



(13) 



These icf 's depend on the ionic and total abundances of oxy- 
gen and s ulfur. A m ap for the total abundance of oxygen was pre- 
sented in iPaper lit We utilized a constant sulfur abundance de- 
termined as the mean of those presented in Paper III. Regarding 
the ionic abundances of oxygen and sulf ur, they were derived 
as part of the work presented in IPaper III However, maps were 
not included there and are displayed in Fig. QT| for complete- 
ness. An extra assumption is needed to calculate the icf's in 
those areas where no measurement for the ionic abundance is 
available. For that, we utilized the information in the other spax- 
els to fit a first-degree polynomial to the relation between the 
[S n],l/l6717,6731/Ha< line ratio and icf(He): 

icf(He) K&S 83 = (0.354 ± 0.012)[Sn]/Ha + (1.033 + 0.002) (14) 

icf(He)p &TP 77 = (0.447 ± 0.026) [S n] /Ha + (1.054 + 0.003) (15) 

icf(He) M13 = (0.763 ± 0.025) [Sn] /Ha + (1.057 + 0.004) (16) 

The maps for the three icf(He) estimations are presented in 
the left column of Fig. [12] In all three cases, the icf is smallest at 
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Fig. 12. Left: Maps for the icf( He) as est i mated from 
the expressions proposed b y Kunth & Sargent (119831 top), 
IPeimbert & Torres-Peimbertl ( 1977 center) and assuming no 
gradient in helium abundance bottom. Right: Corresponding he- 
lium abundance maps. 



the peak of emission for the ionized gas and increases outwards, 
following the ionization structure. 

The icf(He) values are close to one another but variable 
within a range ~ 1.04 - 1.09, ~ 1.06-1.12 and ~ 1.09 - 1.20 
when using Eqs.[TT|[T2] andQ~3] respectively. The corresponding 
maps with the total helium abundance are presented in the right 
column of this figure, while the dependence o n the excitation is 
resented i n Fig . [I0l Note that, when using the Kunth & Sargent 
19831) and IPeimbert & Torres-Peimbertl dl9~77l) icf's. the slope 
of the relation between y and [Oiii]/t5007/Hy6 is slightly posi- 
tive within the errors of the fit. However, according to the crite- 
rion presented for the y + vs. [O ra]/l5007/Hy6 relation, the derived 
helium abundance for these icf's would be effectively consis- 
tent with no positive gradient. The mean total helium abundance 
ranges between 0.086 and 0.091, depending on the assumed icf. 
This range is larger by a factor ~2 than the uncertainties due to 
the measurements (assumed to be traced by the standard devia- 
tion). This implies that the main source of uncertainty is still in 
the assumptions taken on the way to the derivation of the final 
abundance and highlight how, even with a level of data quality 
as high as those utilized here (i.e. data allowing to determine 
locally the physical conditions of the gas, the fluxes of several 
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helium lines and estimations of the contribution of underlying 
stellar populations), achieving an uncertainty <1% is extremely 
difficult. 

The main conclusion of Fig. Qj] is that the relation between 
the total helium abundance, y tot , and the excitation of the gas, 
as traced by [Om]/I5007/Hy8 is consistent with a lack of gra- 
dient in helium abundance. However, if the extra helium were 
produced by the Wolf-Rayet stars, the required amount to be 
detected would be tiny in comparison to the pre-existing he- 
lium. This implies that a positive slope in Fig. Qj] consistent 
with this enrichment would be indistinguishable of the presented 
fits. Exploiting these data to their limits, we can derive the av- 
erage excess in N/He, by comparing the mean abundances in 
the N-enriched and non N-enriched areas. A map for the N/H 
abundanc e can be derived using the O/H and N/O maps pre- 
sented in Paper II. Assuming a value of log(N/0 = -1.3) as 
the limit above which the interstellar medium is enriched in ni- 
trogen, we find a mean N/H of 8.4 x 10~ 6 and 17.4 x 10~ 6 in 



the spaxels without and with enrichment. This implies a flux 
weighted excess in nitrogen of N/H exc = 8.9 x 10~ 6 . Proceeding 
in the same manner with the different estimations of the total 
helium abundance, we find mean He/H values ranging between 
8.49 x 10~ 2 and 9.05 x 10~ 2 for the non-enriched spaxels and 
between 8.68 x 10~ 2 and 9.14 x 10~ 2 , depending on the as- 
sumed icf(He). Differences in helium abundances between the 
enriched and non-enriched zones are ~ 0.9 - 1.9 x 10 -3 . The 
lowest value was derived using our largest icf(He) (i.e. Eq. [TBI 
which, by construction, was defined to minimize any helium 
abundance gradient. The largest value compares well with the 
stella r atmospheric N/He ratios for an N-type Wolf-Rayet star 
fe.g. lSmiffi&WillisLll982h . However, it is also at the limit of the 
uncertainties (~ 2 x 10~ 3 assumed to be traced by the standard 
deviation for the spaxels under consideration in each separate 
group). Therefore, these data appear to be marginally in accord 
with the hypothesis of a putative enrichment of helium due to 
the Wolf-Rayet population in the main H n region, but more im- 
portantly, they stress the difficulties in pushing this methodology 
further in order to confirm this contamination without doubt. 



4. Discussion 

4.1. Kinematics and radiative transfer effects 

The standard approach in the literature to estimate radiative 
transfer effects assumes a negligible influence of the movements 
of the gas in the Hn region, i.e . a static nebula (u> = v/v t h = 0, 
in the formalism presented by iRobbinsl I1968I) . This is the ap- 
proach utilized in Sec. I3.2I The general use of this assumption 
is motivated by the difficulty of obtaining data of sufficient qual- 
ity to allow, on top of measuring the flux of several helium lines 
with good signal to noise, tracing their profiles and clearly iden- 
tify the different kinematic components in a consistent manner 
in all of them. However, it is not unusual for starburst galax- 
ies or H ii regions to present velocity gradients and/or relatively 
high velocity dispersion that can be attributed to outflows (or ex- 
pa nding st ruct ures). In p articular, the kinematic study presented 
in [Paper I and Paper III showed that movements in this H n re- 
gion are significant and can indeed be attributed to an outflow 
ca used by the two embedded Super Star Clusters. Specifically, 
in lPaper II we detected a relatively static (i.e. small velocity gra- 
dient) narrow component on top of a much broader component 
(cr ~ 20-25 km s~') with a velocity gradient of Av ~ 70 km s . 
An additional compo nent was also detected with the higher res- 
olution GMOS data (Pa per II j) . Therefore, these data provide a 
very good opportunity to explore what is the impact of the kine- 
matics in the derivation of the optical depths from an empirical 
point of view. 

Since this is only an exploratory analysis and a consistent 
approach based on the identification and modeling of multiple 
kinematic components in several emission lines would be more 
complex (and not supported by the signal-to-noise of all the He i 
data) than the one presented in Sec. I3.2I we opt for a simpler 
analysis based on the He 1/16678 and He 1/17065 lines. These were 
the pair of lines utilized to estimate t(3889) in the red set. Also 
both were observed with the same instrumental set-up, which 
was the same as for Ha in the multicomponent analysis pre- 
sented in lPaper j . Therefore, we could link on a spaxel- by-spaxel 
basis the multicomponent analysis utilized in iPaper II to simi- 
lar components in He 1/16678 and He 1/17065. The contribution to 
the total flux from each component varies from spaxel to spaxel 
as well as with the line in consideration. Specifically, the broad 
component presents a larger contribution (i.e. ~ 35 - 60%) to 
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Fig. 14. r(3889) f or the tw o main kinematic components pre- 
sented in Fig. 19 oflPaoer il derived from the fits to relations pre- 
sented by Ro bbie (fl968T) for 7,,= 10 000 K. Left: Narrow kine- 
matic component and o» = 0. Right: Broad kinematic component 
and a) — 3. The r(3889) map for this component using the rela- 
tion for a) — (not shown) displays the same structure but with 
values ~ 4 - 6 times smaller. 



He 1/17065 than to He 1^6678 (~ 30 - 50%). This implies differ- 
ent correction factors, / T (7065), for the narrow and broad com- 
ponent that can be converted into optical depths as in Sec. 13.21 
Given the velocity gradient and the velocity dispersion measured 
for the broad component, values for a> — v/v t h — 3 are more ap- 
propriate for this component. Fitting the s ame func ti onal form 
as in Eqs. IBTTOl to the values reported in Rob binsl (1 1968b for 
o) = v/\>th = 3 and T e = 10000 K, we derived the following 
relation: 



/ T (7065L =3 = 1 + 0.279V 



.0.521 



(17) 



The structure of the derived optical depth for both the nar- 
row and the broad components is shown in Fig. [14] and several 
conclusions can be drawn from this figure. Firstly, a comparison 
between the left panel of this figure with those presented in Fig. [8] 
highlights the similarity between the optical depth maps derived 
for the narrow component and for the line integrated analysis. 
To our knowledge, this is the first empirical evidence support- 
ing the traditional assumption that movements in extragalactic 
H ii regions with a known outflow have a negligible effect in the 
estimation of the global optical depth. 

Secondly, a comparison between the two maps in Fig. [14] 
shows how the structure of the two t(3889) is different: while 
the narrow component t(3889) displays a clear peak associated 
with the peak of emission for the GH iiR and decreases outwards, 
the broad component remains relatively high over a large area of 
~ 4" x 2" (~74 pcx37 pc) centered on the double cluster. 

Thirdly, the comparison between the specific values dis- 
played in these two graphics shows that the broad component 
suffers from radiative transfer effects to a larger degree than the 
narrow one. This is a consequence of the adopted relation be- 
tween / T (7065) and t(3889), which is motivated by our kine- 
matic results. If we had used Eq.[5](i.e. assuming cj — for the 
broad component), we would have obtained comparable values 
for the t(3889) in both components (although, still, with a dif- 
ferent structure since this is determined by the f T (1065) itself, 
which in turns depends on the relative flux between He 1/16678 
and He 1/17065). 

There is certainly room for improvement in this analysis. 
For example, although the multicomponent analysis treats in a 
consistent way the line profiles of He 1/16678 and He 1/17065, the 
same physical conditions (i.e. T e , n e , extinction) had to be as- 
sumed for both components. Moreover, we used here only two 



helium lines, in contrast with the more canonical methodology 
for determining helium abundances based in a larger set of he- 
lium lines. These improvements suggest a new approach to the 
study of the helium content in galaxies with IFS. Specifically, 
the 2D analysis presented in the previous sections together with 
the separate analysis of different kinematic components along 
the line of sight (i.e. in a given spaxel) constitute a 3D view of 
the radiative transfer effects in the nebula. 

4.2. Relation between collisional and radiative transfer 
effects and other properties of the ionized gas 

We have evaluated locally both collisional and radiative transfer 
effects in the helium lines, as well as several important physi- 
cal conditions (n e , T e , excitation) and chemical properties (i.e. 
relative abundance N/O) of the ionized gas. Therefore, the cen- 
tral part of NGC 5253 constitutes a good case study to explore 
whether there is a (strong) spatial relation of collisional and/or 
radiative transfer effects with other gas properties. 

As discussed in Sec. 13. II theory predicts a strong dependence 
of the collisional effects on the electron density and, to a lesser 
extent, on the electron temperature. This implies a strong spa- 
tial correlation with the these properties. Indeed, when all the 
spaxels in the area under study are considered, the relations be- 
tween C/R(A1065) vs. n e and T e have Pearson's correlation co- 
efficients of 0.99 and 0.71, respectively. Regarding the other two 
properties, Fig. [I5]presents the relation between the collisional 
effects (as traced by C//?(7065)) and the excitation (as traced by 
[Oiii]/15007/Hy6) and the relative abundance of nitrogen: N/O. 
The excitation presents a degree of correlation similar to n e . The 
correlation for N/O is not as strong as that for n e but still com- 
parable to that for T e . 
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Fig. 15. Relation between collisional effects as traced by the C/R 
ratio for the /17065 line and the excitation {left) and relative abun- 
dance of nitrogen, N/O (right). The Pearson's correlation coeffi- 
cients are indicated in the corners of the individual plots. 

Likewise, Fig. Q~6] shows the relation between radiative ef- 
fects, as traced by r(3889) and different properties of the gas. 
As expected, there is a good correlation with electron density. 
For electron temperatures lower than ~ 10 700 K, Hei emission 
lines do not suffer from radiative transfer effects in a significant 
way. However, above this threshold there is a clear correlation 
between the temperature and the relevance of the radiative trans- 
fer effects. 

More interestingly, as happened with the collisional effects, 
there is a good correlation between the excitation and the rele- 
vance of radiative transfer effects (lower left panel in Fig. HBl . 
It is not clear why the excitation should correlate with the con- 
tribution of collisional effects (mainly) or the radiative transfer 
effects (to a lesser extent). 

Something similar occurs with the relative abundance, N/O. 
At the typical \og(N/0) values for this galaxy (< - 1.3), Hei ra- 
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Fig. 16. Relation between the self-absorption effects as traced 
by t(3889) and different physical and chemical properties in 
NGC 5253. From left to right and from top to bottom, these are 
electron density and temperature, excitation and relative abun- 
dance of nitrogen, N/O. The correlation coefficients are indicated 
in the corners of the individual plots. 

diative transfer effects appear to be negligible. However, those 
locations displaying extra nitrogen appear to show significant 
radiative transfer effects. Rather than a discernible correlation 
between individual spaxels, there appears a correlation between 
the upper limit to the contribution of the radiative transfer effects 
and the amount of extra nitrogen. To our knowledge, there is no 
reason to link these two properties from the theoretical point of 
view. The fact that both the relations with T e and N/ O present 
bimodal behavior (i.e. with "turning points"), and that there is a 
similar structur e for the electron temperature and relative abun- 
dance (Pa per Hi) , implying a good correlation (i.e. r — 0.91) be- 
tween T e and N/ O, points towards the local electron temperature 
as a common cause. 

5. Conclusions 

This is the fourth in a series of articles that make use of IFS- 
based data to study in detail the 2D physical and chemical 
properties of the gas in the main GHuR of the nearby BCD 
NGC 5253. The main goal of this article was to estimate the 
contribution of the collisional and radiative transfer effects on 
the helium emission lines and to map, in turn, the helium abun- 
dance. The major conclusions can be summarized as follows: 

1. The collisional effects on the different helium transitions 
have been mapped for the first time in an extragalactic ob- 
ject. As expected, they reproduce the electron density struc- 
ture. They are negligible (i.e. ~0. 1-0.6%) for transitions in 
the singlet cascade while relatively important for those tran- 
sitions in the triplet cascade. In particular, they can contribute 
up to 20% of the flux in the He 1/17065 line. 

2. The contribution of the collisional effects is sensitive to the 
assumed T e for helium. Specifically, we found differences 
~ 25 - 35% for the ,13889 and A1065 lines two reasonable 
assumptions for the T e sensitivity. Relative differences for 
the other lines were larger. However, this does not have im- 
portant consequences since the contribution of the collisional 
effect to the observed spectrum for these lines is negligible. 

3. We present a map for the optical depth at /13889 in the main 
GHuR of NGC 5253. r(3889) is elevated over an extended 



and circular area of ~30 pc in diameter, centered at the Super 
Star Cluster(s), where it reaches its maximum. 

4. The singly ionized helium abundance, y + , has been mapped 
using extinction corrected fluxes of six He i lines, realistic 
assumptions for T e , n e , and the stellar absorption equivalent 
width as well as the most recent emissivities. We found a 
mean(± standard deviation of 10 3 y + ~ 80.3(+2.7) over the 
mapped area. 

5. We derived total helium abundance maps using three pos- 
sible icf(He)'s. The relation between the excitation and the 
total helium abundances is consistent with no abundance gra- 
dient. Differences between the derived total abundances ac- 
cording to the three methods are larger than statistical errors 
associated with the data themselves, emphasizing how un- 
certainties in the derivation of helium abundances are domi- 
nated by the adopted assumptions. 

6. We illustrated the difficulty of detecting a putative helium en- 
richment due to the presence of Wolf-Rayet stars in the main 
GH nR. This is due to the comparatively large amount of pre- 
existing helium. The data are marginally consistent with an 
excess in the N/He ratio in the nitrogen enriched area of the 
order of the atmospheric N/He ratios in W-R stars. However, 
this excess is also of the same order of the uncertainty esti- 
mated for the N/He ratios in the nitrogen enriched and non- 
enriched areas. 

7. We explored the influence of the kinematics in the evalua- 
tion of the He i radiative transfer effects. Our data empirically 
support the use of the traditional assumption that motions in 
an extragalactic H n region have a negligible effect in the es- 
timation of the global optical depths. However, individually, 
the broad kinematic component (associated with an outflow) 
is affected by radiative transfer effects in a much more sig- 
nificant way than the narrow one. 

8. The local relationships between the contribution of colli- 
sional and radiative transfer effects to the helium lines and 
different physical and chemical properties of the gas have 
been explored. Interestingly, we found a relation between the 
amount of extra nitrogen and the upper limit of the contribu- 
tion from radiative transfer effects that requires further in- 
vestigation. We suggest the electron temperature as perhaps 
a common agent causing this relation. 
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